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PREFACE 

With the formal exercise in solving the types of ordinary 
differential equations that usually occur it is the object of 
this text to combine a thorough drill in the solution of prob- 
lems in which the student sets up and integrates his own 
differential equation. For this purpose certain topics in 
mechanics and physics needed in groups of problems are 
briefly presented in the text. 

The problems have been collected from a variety of sources 
among which the author wishes particularly to mention the 
Advanced Calculus of Professor E. B. Wilson and the notes 
on Mathematics for Chemists prepared by Professors W. 
K. Lewis and F. L. Hitchcock. 

H. B. PHILLIPS 

Cambridge, Mass., 
Feb. 15, 1922. 
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DIFFERENTIAL EQUATIONS 

CHAPTER I 

FIRST ORDER EQUATIONS, VARIABLES 

SEPARABLE 

1. Definitions. — In this chapter we consider problems in- 
volving two variables one of which is a function of the other. 
It is often possible from the statement of a problem to ob- 
tain an equation involving the difiFerentials or the derivatives 
of the variables. Such an equation is called a differential 
equation. Thus 

{p^ + 2/^) dx + 2xydy ^ 
and 

dx^ dx 

are differential equations. 

A solution of a differential equation is an equation con- 
necting the variables such that if the derivatives are cal- 
culated from it and substituted in the differential equation, 
the latter will be satisfied. Thus 

y = x^ — 2x 

is a solution of the second equation above; for when x^ — 2x 
is substituted for y the equation is satisfied. 

An equation containing only first derivatives or differ- 
entials is called an equation of the first order. In general, 
the order of a differential equation is the order of the highest 
derivative occurring in it. 

2. Separation of the Variables. — If a differential equa- 
tion has the form 

fi(x)dx+f2(y)dy = 0, (2a) 

1 
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2 DIFFERENTIAL EQUATIONS Chap. I 

one term containing only x and dx, the other only y and dt/, 
the variables are said to be separated. The solution is ob- 
tained by integration in the form 

ffi {x) dx + Ju (y) dy = c, (2b) 

• 

where c is a constant of integration. 

Since the integration formulas contain a single variable, 
if the variables are not separated, we cannot solve the 
equation in this way. Thus, if 

xdy + (1 " y) dx = 0, 

since x dy cannot be integrated, we cannot obtain a solution 
by direct integration. By division we can however reduce 
this equation to the form 

^^ +^=0 (2c) 



1-2/ X 

in which the variables are separated. The solution is then 

In X — In (1 — 2/) = c. 

When the variables can thus be separated the differential 
equation is called separable. An equation of the form 

Mdx + Ndy = 

is separable when each of the coefficients Af , iV is a function 
of only one variable or the product of factors each contain- 
ing a single variable. 
3. Different Forms of Solution. — The solution 

' ^l. • hi X - hi (1 - 2/) = c (3a) 



can be written 



whence 



ln= = c, 

1-2/ 



^ = e* = fc. 

1-2/ 
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Since c is an arbitrary constant, k is also arbitrary. The 
solution could then be written 

x^c(l-y) (3b) 

where c is an arbitrary constant. It could also be written 

1 " y — ex (3c) 

y -1 = ex. (3d) 

Any one of the equations (3a), (3b), (3c), (3d) is the solution 
of (2c), but of course the constant has a different meaning 
in each case and so two of these could not be used simultane- 
ously. 

Example. Solve the equation {1 + x^)dy — xy dx = 0. 
Separating the variables, this becomes 

dy xdx _ ^ 
2/j 1 + a^ ~ "' 

whence 

In y — i In (1 + x^) = const. 

Since any constant is the logarithm of another constant, 
this can be written 

In 2/ — i In (1 + a:^) = In c, 
whence 

y = cVl + x\ 

This answer. could equally well be written in any one of the 
forms 

2/2 = c2 (1 + x2), 

2/2 = c (1 + x% 
ey^ = 1 + a?. 

EXERCISES 

Solve the followmg equations: 

1. tan X sin* ydx + cos* x cot ydy = Q. 

a. (a^ + x)dx + {y — x*2/) dy = 0. 
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8. (a^* + x)dx + (a^ — y) dy = 0. 

4. Derivative Relations. — In many cases one or more 
of the quantities occurring in a problem is a derivative. 
An equation satisfied by these quantities is then an equation 
containing a derivative, i.e., a differential equation. 

Thus it may be known that the slope of a curve is a given 

function of x and y. Since the slop)e is -^ > the curve can be 
obtained by solving the differential equation 

g = /(x,y). 

Again, it may be known that the velocity oi a moving 
particle is a given function of the distance s and time t. 
The differential equation is then 

More generally, if the rate of change of a quantity z is 
known to be a function / (x, t)j then 

6. Determination of Constants. — Since the constant of 
integration may have any value whatever, there are an 
infinite niunber of solutions of a given differential equation. 
A pair of corresponding values of the variables is however 
usually known. By substituting these in the solution the 
constant can be determined and so a definite solution be ob- 
tained. 

In many cases the derivative is known merely to be pro- 
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portional to a ceptain function / (x, y). The differential 
equation is then 

where k is constant. If two pairs of corresponding values 
Xi, y\ and X2, 2/2 are known, by substituting in the solution 
both k and c can be determined. 

The statement of a problem thus consists of two parts. 
One part contains conditions true at all places or times. 
From this the differential equa- 
tion is determined. The second 
part contains conditions true at 
a single place or time. These 
are used to determine the con- 
stants. 

Example 1. Find the curve 
passing through (2, 3) such that 
the part of the tangent between 
the coordinate axes is bisected 
at the point of tangency. 

Every tangent is bisected at the point of tangency. Let 
P (x, y) be the middle point of the tangent AB. Then by 
similar triangles . ,, ^^^^ ^ e». 

OA = 2t/, OB = 2x. ' / ^ . * 




Fig. 5. 



The slope of the curve at P {x, y) is 



dy ^ _0A ^ __y 
dx OB x' 



\i 



<r. \ ^ '^ 



1^ 



This can be written 



I ' «- 'X' 






T , ; 'x t ■• V 



X 



• < ' ' 



the solution of which is 



xy = c. 
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Since the curve passes through (2, 3), we must then have 

2 (3) = c. 

Hence the equation of the ciu^e is 

xy= 6. 

Example 2. Radiiim decomposes at a rate proportional 
to the amount present. If half the original quantity disap- 
pears in 1600 years, what percentage disappears in 100 years? 

Let R be the amount of radium present at time t. The 

dR 
rate of decomposition is measured by — -^ • Since this is 

proportional to R^ 

dt *^' 
where k is constant. Hence 

and 

hi R = kt + c. 

Let Ro be the amount at the start. Substituting ^ = 0, 
jB = jBo, we have 

In Ro = c. 

Substituting this value of c and transposing, we have 

In^ =kt. 

When t = 1600, 72 = § JBo. Hence 

hi J = 1600 fc, 
whence 

1600 
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When t = 100 we therefore have 



which gives 



^|=-S-l«^=--«*33 



'^~ aVIdO* 



This shows that 95.8% remains at the end of 100 years 
and so 4.2 % disappears. 

6. Differential Relations. — It is usually easier to find 
relations between the first differentials of the variables than 
between the variables themselves. This is due to certain 
simplifying assumptions that may be made without affecting 
the results. Thus, so far as first differentials are concerned, 
a small part of a curve near a point may be considered 
straight and a part of a surface plane; during a short time 
dt a particle may be considered as moving with constant 
velocity and any physical process as occurring at a con- 
stant rate. The reason these assumptions give a correct 
result is because the ratio of differentials is by definition 
the limit of the ratio of increments, and as the increments 
approach zero these simple conditions become more and 
more approximately satisfied. 

Methods of setting up differential relations in this ap- 
proximate way are often called differential methods. As 
here stated these methods apply only to first differentials, 
or first derivatives. A correct equation containing second 
derivatives would not in general be obtained, by considering 
a small part of a curve as straight and in a small interval 
a physical process as occurring at a constant rate. 

Example 1. Find the shape of a mirror such that all light 
coming from one fixed point is reflected to another fixed 
point. 

Let the light from F (Fig. 6) be reflected to F\ The 
mirror must have the form of a surface of revolution. Other- 
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wise light passing out from F in a, plane through FF' would 
not be reflected in the same plane and so could not go to F\ 
Let PQ be an infinitesimal arc. With F as center construct 

the arc QS and with F' 
as center, the arc QRi 
Consider PQS and PQR 
as right triangles. They 
have the common hypo- 
tenuse PQ. Also, since the 
angles of incidence and 
reflection are equal, 

Z QPS = z QPR. 




Fig. 6. 



Hence the right triangles are equal and so 

PS = PR. 



(6a) 



Let r = FP, r' = F'P. In passing from P to Q the in- 
crease in r is 

dr = PS' 

and that in r' is 

dr' = - PR. 

Hence, from (6a), 

dr = — dr' 

and so 

r + r' = const. (6b) 

The section of the mirror by any plane* through FF' is there- 
fore an ellipse with F, F' as foci. 

Example 2. The sum of $100 is put at interest at 5% per 
anniun under the condition that the interest shall be com- 
poimded at each instant. How many years will be re- 
quired for the amount to reach $200? 

Let A be the amoimt at the end of t years. In the short 
time dt the increase will equal the interest 

(U = .05 A di. 
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Integrating between the limits A = 100 and A = 200, 
we get 

^ = .05« 
100 A 

whence 



/ 

«/ia 



, 1 , 200 ,-_ 

''=:65'°ioo = ^^-^y^- 

7. Flow of Water from an Orifice. — If there were no 
loss of energy, the velocity with which water would issue 
from an orifice at depth h below the surface would be that 
acquired by a body in faUing the distance ft, namely, 

. . V2gh. 

Because of friction and the converging form of the stream, 
the average velocity with which the water actually issues is 

V = c y/2 gh, 

where, for ordinary small orifices with sharp edges, 

c = 0.6 

approximately. 

Example. Find the time 
required for a hemispherical 
bowl 2 ft. in diameter to empty 
through an inch hole at the FiiTV 

bottom. 

Let h be the depth of water at time t and 

r = Vl - (1 - hy 

the radius of the circle forming its surface. The water 
which issues in time dt generates a cyUnder of altitude v dt 
with an inch circle as base. Its volume is 




^(24) ^"^^^ 
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This causes the loss from the surface of a slice of water with 
radius r, altitude — dh, and volume 

-7rr^dh=- t Qv" -2h)dh. 
Hence 

T (A* -2h)dh = 'Q)*"'" = '^1^|^<* 

Separating the variables, we have 

t^ 120 r (At - 2 A») dh = 112 sec. 

8. Equation of Continuity. — In a physical process there 
may occur an element which is neither created nor destroyed. 
The amount of this element in a given region only changes 
when some comes in or goes out through the boundary. 
In such a case the obvious equation 

increase = income — output 

is sometimes called an equation of continuity. Stated in 
diflferential form this may give a differential equation by 
which the variation of the particular element can be de- 
tennined. 

The concentration c of a particular substance is the amount 
of that substance in imit volimie. If the concentration is 
uniform the amount in the volume v is then cv. 

Example, In a tank are 100 gallons of brine containing 
60 lbs. of dissolved salt. Water runs into the tank at the 
rate of 3 gals, per minute and the mixture rims out at the ^ 
same rate, the concentration being kept uniform by stir- 
ring. How much salt is in the tank at the end of one hour? 

Let X be the amount of salt in the tank at the end of t 
minutes. The concentration is then 

c = jTT^ pounds per gallon. 
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In the time dt, 3 dt gals, of water come in and 3 dt gals, of 
brine containing 3cdt lbs. of salt go out. Hence the change 
in the amount of salt in the tank id 

dx = —Scdt = — r^ dt. 

The amount of salt at the end of one hour is then deter- 
mined by 

rdx^ 3_ p , 

J^x lOOjo ' 

whence 

and 

X = 8.27 lbs. 

9. Flow of Heat. — If the temp)eratures at the boimding 
surfaces of a body are kept constant, the body will ulti- 
mately approach a steady state in which the temperatures 
at different points may be different but the temperature at 
a given point no longer changes with the time. In many 
cases the temperature T is a function of a single coordinate 
X. By Newton's law, the rate at which heat flows across an 
area A perpendicular to x is then 

-fcA^=(3, (9a) 

where fc is a constant called the conductivity of the material. 
If we have a series of surfaces A such that the heat flow- 
ing across one flows across all the others, the equation of 
continuity has the form 

Q = const. (9b) 

If then A is expressed in terms of x, the solution of (9a) 
gives 7 as a function of x. By substituting the values of 
X and T at two boundaries, the constant of integration and 
Q can then be determined. 
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Example. A hollow spherical shell of inner radius 6 cm. 
and outer radius 10 cm. is made of iron (k = 0.14). If the 
inner temperature is 200° C. and the outer 20° C, find the 

temperature at distance r from 
the center and the amount of 
heat per second that flows 
outward through the shell. 

By symmetry the flow . of 
heat is seen to be radial. At 
distance r from the center the 
area across which heat is flowing 
is the spherical surface 

A = 4 7rr2. 
Fig. 9 

Since there is no accumulation 
of heat between the surfaces, the same amount flows across 
ea€h spherical surface and so equation (9) is 

— 4:TrkT^^- = Q = const. 
dr 

Separating the variables and integrating, we get 

4TkT = ^ + c. 
r 

Substituting T = 20, r = 10 and T = 200, r = 6, we find 

C = - 1000 Trfc, Q == 10,800 rk, 

r=?^-250. 



The rate of flow through the shell is 

Q = 10,800 wk = 4750 cal./sec. 

10. Second Order Processes. — In a problem con- 
taining two independent variables x and y it is sometimes 
stated that a quantity z is proportional to x and also pro- 
portional to y. What is meant is that, y being constant 
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z is proportional to x, and x being constant z is proportional 
to y. Both statements are expressed by the equation 

Thus, tne rate at which a substance x dissolves is pro- 
portional to the amount of x present. It is also propor- 
tional to the difiFerence between c, the concentration of x in 
the solvent, and s its concentration in a saturated solution. 
These statements are both expressed by 

T7 = fcx (s — c). (10a) 

The fact that the total amount of x present (solid and in 

solution) is constapt gives an equation of continuity from 

SiX 
which we can express c in terms of x and so express -jj as a 

quadratic function of x, A process in which the rate of 
change of x is a quadratic fimction of x is called a second 
order process. 

Example. Sulphur is to be removed from an inert 
material by extraction with benzol. By using a large 
amount of benzol it is found that half the sulphur can be 
extracted in 42 min. If the material contains 6 gms. of 
sulphur and 100 gms. of benzol are used, which if saturated 
would dissolve 11 gms. of sulphur, how much of the sulphur 
will be removed in 6 hrs.? 

Let X be the amoimt of sulphur undissolved at time U 
The concentration of sulphur in a saturated solution is 

^ = ^gms. sulphur per gm. benzol. 

The differential equation for x is then 

%=kx (0.11 - c). (10b) 
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If we use a very large amount of benzol, c will be very smaU 
and so this may be replaced by 

I =*^ (0.11). 

Since x then varies from 6 to 3 in 42 min., 

r—^OAlkT'dt, 

•76 X Jo 

whence 

fc = - 0.16. 

If now we use 100 gms. benzol, when there are x gnis» 
sulphiu" imdissolved, there will be 6 — a; in solution and so 

6 - X 

"="Ior- 

Hence (10b) can be written 

~ 0.15 a:(o.ll - ^j^) = - .0015 x (5 + x), 

and 

dx r^ 



r dx _ _ ^^^p 
Ja x(x+ 5) ' Jo 



dt 



1 11a; _ ^^ 
^''6(x + 5)" "^''^ \ 

which gives x = .19 gms. as the amount of sulphur undis- 
solved at the end of 6 hrs. 

PROBLEMS 

1. Find the equation of the curve passing through the origin such 
that the part of every tangent between the avaxis and the point of 
tangency is bisected by the y-&ida. 

2. Find the curve passing through the point (2, 0) .such that the 
part of the tangent between the y-axis and point of tangency is of 
length 2. 

8. If in the culture of yeast the amount of active ferment doubles 
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in one hour, how much may be anticipated at the end of 2} hours at 
the same rate of growth. 
/ 4. If the activity of a radioactive deposit is proportional to its 
rate of diminution and is found to decrease to } its initial value in 4 ^ 
days, find the value of the activity as a function of the time. 

6. The retarding effect of fluid friction on a rotating disk is pro- 
portional to its angular velocity. If the disk starts with a velocity of 
100 revolutions per minute and revolves 60 times during the first 
minute, find its velocity as a function of the time. ^ 

6 6. According to Newton's law the rate at which a substance cools 
in air is proportional to the difference of its temperature and that of i 
the air. If the temperature of the air is 20° C. and the substance cools 
from lQp° to 60° in 20 min., when will its temperature become 30°? 

7. A substance is undergoing transformation into another at a rate 
proportional to the amount of the substance remaining untransformed. 
If that amount is 31.4 at the end of 1 hr. and 9.7 at the end of 3 hrs., 
find the amount at the start and find how many hours will elapse 
before only 1% will remain. 

8. When a liquid rotates about a vertical axis, show that its smface 
forms a paraboloid of revolution. Observe that the weight of a par- 
ticle of water at the smface and its centrifugal force must have as 
resultant a force perpendicular to the surface. 

9. Through each point of a curve lines are drawn parallel to the axes 
to form a rectangle two of whose sides lie on the axes. Find the curve 
which cuts every rectangle of this kind into two areas one of which is 
twice the other. 

10. Find the siuface of a mirror such that all light from a fixed 
point is reflected parallel to a fixed Une. Take the fixed point as or- 
igin 'and let the light be reflected parallel to the a;-axis. Use the polar 
coordinate r and the x-coordinate as variables. 

11. In a certain type of reflecting telescope light converging toward 
a fixed point is reflected by a mirror to another point. Find the shape 
of the mirror. 

12. If a man can earn 5 dollars per day over expenses and keep 
his earnings continuously invested at 6% compound interest, how 
long will it take to save $25,000? 

13. If a man can earn a dollars per year over expenses and keep his 
savings continuously invested at 6% compoimd interest, how long will 
be needed to obtain an income of a dollars per year from investments? 

- 14. The amount an elastic string of natural length { stretches under 
a force F is WF, k being constant. - Find the amount it stretches when 
suspended from one end and allowed to stretch imder its own weight w. 
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16. Find the amount the string of the preceding problem stretches 
if it is hung up with a weight P attached to the lower end. 

16. Consider a vertical column of a\r and assume that the pressure 

^ at any level is due to the weight of air above. Find the pressure as a 

function of the height if the pressure at sea level is 14.7 lbs. persq. 

in., and at an elevation of 1600 ft. is 13.8 lbs. per sq. in. Assume 

B^jje's law that the density of the gas is proportional to the pressure. 

ll.j If air in moving from one level to another expands without 
receiving or giving out heat, that is, adiabatically, 

where p is the pressure, p the density, and k, n constants. Assuming adia- 
batical expansion, if the density at sea level is .08 lbs. per cu. ft. and 
the pressure 2100 lbs. per sq. ft., find the height of the atmosphere. 
18. If the coeflficient of friction between a belt and pulley is /a and 
the angle of lapping a, show that the tensions 7'i, Tt in the two sides 
of the belt when it is slipping satisfy the equation 

fl9. If the velocity is high .centrifugal force reduces the pressure of 
the belt upon the pulley. A^uming the weight of the belt to be t& lbs. 
per unit length and its velocity t;, find the equation connecting the 
tensions in the two sides of the belts. 

20. The end of a vertical shaft of radius a is supported by a flat- 
step bearing. If the horizontal surface of the bearing carries a 
uniform load of p lbs per sq. in. (new bearing) and the coefficient of 
friction is At, find the work done against friction in one revolution. 

An old bearing is worn a little more at the edge than at the center. 
Ultimately the pressure varies in such a way that the wear is the same 
at all points. Assmning that the wear at any point is proportional to 
the work of friction per unit area at that point, find the law of vari- 
ation of pressure and show that with a given total load the work per 
revolution is only J that in a new bearing. 

21. Assuming that the density of sea water imder a pressure of p 
lbs. per sq. in. is 

1 + 0.000003 p 

times its density at the surface, show that the surface of an ocean 5 
miles deep is about 450 ft. lower than it would be if water were incom- 
pressible. 

22. A cylindrical tank with vertical axis is 6 ft. deep and 4 ft. in 
diameter. If the tank is full of water, find the time required to empty 
through a 2-inch hole at the bottom. 
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83. Find the time of emptying if the axis of the tank in the preceding 
problem is horizontal. 

24. Two vertical tanks each 4 ft. deep and 4 ft. in diameter are 
connected by a short 2-inch pipe at the bottom. If one of the tanks is 
full and the other empty, find the time required to reach the same 
level in both. Assume that the velocity through the pipe is the same 
as that through an orifice under the same effective pressure. 

86. Into a tank of square cross-section, 4 ft. deep and 6 ft. in diam- 
eter water flows at the rate of 10 cu. ft. per minute. Find the time 
required to fill the tank if at the same time the water leaks out through 
an inch hole at the bottom. 

26. If half the water runs out of a conical funnel in 2 min., find the 
time required to empty. 

27. A vertical tank has a sUght leak at the bottom. Assuming 
that the water escapes at a rate proportional to the pressure and that 
^ of it escapes the first day, find the time required to half empty. 

28. In a tank are 100 gals, of brine containing 50 lbs. of dissolved 
salt. Water runs into the tank at the rate of 3 gals, per min., and the 
mixture runs out at the rate of 2 gals, per min., the concentration being 
kept uniform by stirring. How much salt is in the tank at the end of 
one hour? 

29. Suppose the bottom of the tank in the preceding problem is 
covered with a mixture of salt and insoluble material. Assume that 
the salt dissolves at a rate proportional to the difference between the 
concentration of the solution and that of a saturated solution (3 lbs. 
salt per gal.) and that if the water were fresh 1 lb. salt would dissolve 
per minute. How much salt will be in the solution at the end of one 
hour? 

30. Oxygen flows -through one tube into a hter flask filled with air 
and the mixture of oxygen and air escapes through another. If the 
action is so slow that the mixture in the flask may be considered uni- 
form, what percentage of oxygen will the flask contain after 5 hters have 
passed through? Assume that air contains 21% oxygen. 

31. The air in a recently used class-room 30' X 30' X 12' tested 
0.12% carbon dioxide. How many cu. ft. air containing 0.04% CO2 
miyst be admitted per minute that 10 minutes later it may test 0.06% 
COt. 

82. If the average person breathes 18 times per minute exhaling each 
time 100 cu. in. containing 4% CO2, find the per cent CO2 in the air of 
a class-room i hour after a class of 50 enters, assuming the air fresh 
at the start and that the ventilators admit 1000 cu. ft. fresh air per 
minute. Let the volume of the room be 10,000 cu. ft. and assume 
that fresh air contains 0.04% COa. 
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A factory 20(y X 45' X 12' receives through the ventilators 
10,000 cu. ft. fresh air per minute containing 0.04% COi. A half hour 
after the help enters at 7 a.m* the COs content has risen to 0.12%. 
'What value is to be anticipated at noon? 

84. A brick wall (k » 0.0015) is 30 cm. thick. If the inner surface 
is at 20° C. and the outer at O^C, find the temperature in the wall as 
a function of the distance from the outer smface. Also find the heat 
loss per day through a square meter. 

36. A steam pipe 20 cm. in diameter is protected with a covering 
10 cm. thick of magnesia (k — 0.00017). If the outer surface is at 
30° C. and the surface of the pipe 160° C, determine the temperature 
in the covering as a function of the distance from the center of the pipe. 
Also determine the heat loss per day through a meter length of the pipe. 

N I 36. A wire whose resistance per cm. length is 0.1 ohm is imbedded 
along the axis of a cylindrical cement tube of radii 0.5 cm. and 1.0 
cm. An electric current of 5 amp. is found to keep a temperature 
difference of 125° C. between the inner and outer surfaces. What is 
the conductivity of the cement? ; 

37. The amount of light absorbed in passing through a thin sheet of 
water is proportional to the amoimt falling on the surface and also 

r proportional to the thickness of the sheet. If one-half the light were 
absorbed in penetrating 10 ft., how much would reach the depth of 
100 ft.? 

38. A porous material dries in a confined space at a rate proportional 
to its moisture content and also to the difference between the moisture 
content of air and that of saturated air. A quantity of material con- 
taining 10 lbs. of moisture was placed in a closed storeroom of volume 
2000 cu. ft. The air at the beginning had a humidity of 25%. Sat- 
urated air at the given temperature contains approximately 0.015 
lbs. moisture per cu. ft. If the material lost half its moisture the first 
day, estimate its condition at the end of the second day. 

39. How long would be needed for the substance of the preceding 
problem to lose 90% of its moisture if the humidity of air is kept at 
25% by ventilation. 

40. A mass of insoluble material contains 30 lbs. of salt in its pores. 
The mass is agitated with 30 gals, of water for 1 hour when one-half the 
salt is found to be dissolved. How much would have dissolved in the 
same time if we had used double the amount of water? Assume the 
rate of solution proportional to the amount of imdissolved salt and also 
proportional to the difference between the concentration of the solu- 
tion and that of a saturated solution (3 lbs. salt per gal.). 

41. A mass of inert material containing 5 lbs. of salt in its pores 
agitated with 10 gals, of water. In 5 minutes 2 lbs. of salt have 
solved. When will the salt be 99% dissolved? 



CHAPTER II 

OTHER FIRST ORDER EQUATIONS 

11. Exact Differential Equations. — The equation 

Mdx + Ndy = (11a) 

is called exact if their exists a function u with total diflPerential 

du = Mdx + Ndy. (lib) 

In this case (11a) can be written 

du = 

and so its solution is 

u = c. 

It is shown in calculus that M dx + N dy is an exact 
differential when and only when 

dM dN ^ ... . 

To determine whether an equation is exact we therefore 
calculate these partial derivatives and observe whether they 
are equal. 

To solve the equation it is necessary to find the function 
u whose differential is M dx + N dy. The terms can often 
be arranged in groups each of which is an exact differential. 
The value of u is then' obtained by integrating these groups 
separately. , 

If this cannot be done, the solution can be determined 
from the fact that 

19 
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By integrating with y constant, we get 

u = I M dx+f (2/), 

the constant of integration being possibly a function of y. 
This function of y can be found by equating the diflferential 
of M to N dx + N dy. Since df (y) gives terms containing 
y only, f (y) can usually be found by integrating the terms 
in N dy that do not contain x. In exceptional cases this 
may not give the correct result. The answer should there- 
fore be tested by differentiation. 

Example 1. (0^ + 2/) dx + (2 y + x) dy = 0. This equa- 
tion can be written 

xdx + 2ydy + (ydx + xdy) = 0. 

It is therefore exact and its solution is 

x^ + 2y^ + fxy = c. 
Example 2. e\ dx + {xf? — 2y) dy = 0, 
In this case 

dy, dy^ 

These derivatives being equal, the equation is exact. Hence 

u=je^dx = x^+f (y)y 

du= ^dx + ^x^ + /'(2/)l dy, 

where f{y) is the derivative of / {y). Comparing with the 
original equation, we see that 

and so 
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The solution is 

x^ — y^ — c. 

The value of / {y) could have been obtained by integrating 
— 2ydy which is the part of 

(xey ''2y)dy 

not containing x. 

12. Integrating Factors. — If an equation of the form 

Mdx + Ndy = 

is not exact it can always be made exact by multiplying by 
a proper factor. Such a multiplier is called an integrating 
factor. 

For example; the equation 

y (1 + xy) dx — xdy = 

is not exact. It can however be written 

y dx — xdy + xy^ dx = 0. 

Dividing by y^, 

ydx — xdy . , ^ 
5 + a; ox = 0. 

y" 

Both terms of this equation are exact differentials. The 
solution is 

An integrating factor in this case is —^ • 

if 

While an equation of the form M dx + N dy ^ Q always 
has integrating factors, there is no general method of finding 
them. 

13. Linear Equations. — A differential equation of the 
form 

%+Py-Q> (13a) 
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where P and Q are functions of x or constants, is called 
linear, A linear equation is thus one of the first degree in 
one of the variables (y in this case) and its derivative. Ajiy 
functions of the other variable may occur. 
When the equation is written in the form (13a); 

Jpdz 
e 

is an integrating factor; for, when multiplied by this factor^ 
the equation becomes 

e % + ye P=e Q. 
The left side is the derivative of 

fPdx 

ye 
Hence 

fPdx pfPdx 

ye = I e Qdx + e (13b) 

is the solution. 

Examvle 1. -^ + -y = of' 

ax X 

In this case 

/Pdx= /-dx = 21na: = lnx*. 

Hence 

fPdx InaJ* 

The integrating factor is therefore a?. Multiplying by 
s? and changing to differentials, the equation becomes 

3i?dy + 2xydx = afidx. 
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The solution is 

x^J/ = gX* + c. 

Example 2. (1 + y^) dx — (xy +y + y^) dy = 0. 
This is an equation of the first degree in x and dx» 
Dividing by (1 + y^) dy, it becomes 

dx y _ 

dy l + y"""^^' 

P is here a function of y and 

fPdy J 

6 = > 

multiplying by this factor, the equation becomes 

dx xydy _ y dy 

whence • 

and 

x^l + y' + c ViT¥' 

14. Equations Reducible to Linear Form. — An equation 
of the form 

% + Py^Qit (14) 

where P and Q are fimctions of x can be made Unear by a 
change of variable. Dividing by.y^it becomes 

2r^ + ^2/-^+' = Q. 
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If we take 

as new variable the equation becomes 



1 <iw , D n 



which is linear. 



Example. g + |y=g. 



Division by y^ gives 



Let 



Then 



whence 



j,dy , 2 « 1 



u = trK 



^dy Idu 



dz 2dx 

Substituting tnese values we get 



and so 



_ldu 2 ^1 
2dx X 7? 



du_4: ^ _2 

dx X 7? 



This is a linear equation with solution 



or, since u = y~^f 



"=0+*^' 



y^ 3 a? 
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16. Homogeneous Equations. — A function / (x, y) is said 
to be a homogeneous equation of the nth degree if 

f(tx,'ty) = P'fix^y). 

Thus, Vx^ + t/^is a homogeneous function of the first 
degree; for 

Vx^fi + y^fi = t Vx* + y\ 

It is easily seen that a polynomial all of whose terms are of 
the nth degree is a homogeneous fimction of the nth degree. 
The differential equation 

Mdx + Ndy = 

* 

is called homogeneous if M and N are homogeneous func- 
tions of the same degree. To solve a homogeneous equation, 
substitute 

or 

X =^ vy. 

The new equation will be separable. 
Example 1. x--^ — y = Vx* + y^. 

This is a homogeneous equation of the first degree. Sub- 
stituting y = vx, it becomes 

x(v + X-^J — vx = Vx^ + t;2^^ 

whence 

This is a separable equation with solution 

x=^ c(v + vr+n^). 
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1/ 
Repladng vhy -f transposing, squaring, etc., 

X 

X* - 2 q^ = c*. 
Example2. y (gj + 2x^ - y = 0. 

Solving for ^> we get 

dy^ -z± Vx^ + y* 
dx y 

whence 

ydy + xdx = ±>/x* + y^ dx. 

This is a homogeneous equation of the first degree. It is 

however much easier to divide by Vx* + y* and integrate 
at once. The result is 

xdx + ydy^^^^ 
Vx* + y2 

Integration gives 

V^a? + y^ = c it X 

and so 

y^ = c* ifc 2 ex. 

Since c may be either positive or negative, the answer is 

equivalent to 

y^ = c^ + 2 ex. 

16. Change of Variable. — We have solved the homo- 
geneous equation by taking 

y 

X 

as new variable. It may be possible to reduce any equa- 
tion to a simpler form by taking some function m of x and 
y as new variable or by taking two functions u and v as new 
nables. 
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If the differential equation only is known some expression 
appearing in the equation may be a good variable. Thus it 
often happens that y appears only in the combinations 

2/*andy^- By taking 

. dy Idu 

a simpler equation is obtained. 

If the equation is obtained in the solution of a problem, 
any quantity which plays a prominent role in the statement 
of the problem may be a good variable. Thus, in solving 
the reflector problem (Art. 6) we used as variables the dis- 
tances from the two points which were directly suggested 
by the problem itself. 

dv 
Example, (x — vY -r- = «*• 

Let X — y = w. Then 

- _dy _ du 
dx dx 

and the differential equation becomes 

vr — or = vTj-' 

dx 

The variables are separable. The solution is 

, a, u — a . 
2 u + a 

whence 

a, x — y — a , 

y = 7: In 2— i he. 

*^ 2 X — y + a 
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17. Simultaneous Equations. — We often have two dif- 
ferential equations 

containing two dependent variables x and y and their de- 
rivatives with respect to the same independent variable t. 
It may be possible to combine the equations algebraically 
so as to get an equation containing only one dependent 
variable, which may be x or y or any function of x and y. We 
solve for this and substitute in one of the original equations 
to complete the solution. 

Example 1. ^ = *i (^ "" V)* ^ ^ ^V* 

The second equation contains only one dependent variable 
y. Its solution is 

y = Cie^. 

Substituting this value in the first equation, 

at 
This is a linear equation with solution 

hi — «i 
Example 2. ^ "= ^ "" 2/> rf? ~ ^ ^* 

Multiplying the first equation by a constant a, the second 
by b and adding, we get 

d 

^(ay + bx) = ox -f- (2 5 — a) y. 

The right side of this equation will be a multiple of the 
expression in parentheses if 

a _ 26 — o 
6 ~ a 



\ 
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A solution of this is 

a = 6 = 1, 

and so 

Solving for x + y, 

x + y = citf. 

Substituting 

X ^ Ci€f -y (17) 

in the first equation we get 

1+2, = «,«•. 

This is a linear equation with solution 

y = Cifi-^ + o<^ief. 

Substituting in (17), we get 

X = -^Cief " oar^. 

EXERCISES 
Solve the following differential equations: 

8. (x» + I/*) dx + 2 xy dy = 0. 

4. («« + I/*) dx - 2 xj/ dj/ = 0. 

6. y dx + (x + J/) dy = 0. 

6. y dx — (x + J/) dj^ « 0. 

^ 7. X dj/ + y dx = ^ dx. 
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10. a^^-2xy-3y. 

n. (2x^ -y)dx+xdy ^0. 

^m. dy 

12. tan ^ ;^ ■" !/ " <»• 

18. (x«- l)l<ij/ + (a^+3xy V5^^ri)(to -0. 

14. yei' (to = (|/» + 2 xe^) dy, 

15. (irj/ 6» + V*) (to — z^dy — 0. 

, 16. g + y-xy. 

17. xg-3y+«*»»-0. 

18. xdx '{•ydy — xdy — ydx, 

19. (a:2/» — x) (to + (y + xy) dy = 0. 

20. xdy — ydx ^ Va:« + y* (to. 

21. (sin X •\-y)dy '\- (ycoBX_— x*) (to = 0. 

22. xdy — ydx ^ X Vx» -i-y^dx, 
28. (1+ x») (ij/ + (a?|^ - x») (to = 0. 

24. y (to = (v* — x) (iy. 

25. y J + y* cot x = cos x, 

26. (x + j/-l)(to + (2»+2y-3)(iy -0. 

27. 3»«g-|/»-x. 

81. (2a;+3y-l)(to + (4x + 6j/-5)(ij/ =a 
(3 1/» + 3 a?|^ + «*) <to = (x» + 2 a?|^) (ij/. 



•« a ^ <iy 1 o • A dx . . 
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PROBLEMS 

1. Using rectangulax coordinates, find the shape of a reflector such 
that light coming from a fixed point is reflected parallel to a fixed line. 

2. In Example 1, Art. 8, suppose the outflow passes through a 
second 100 gal. tank initially filled with pure water. How much salt 
will this tank contain at the end of 1 hour? 

If i is the current, the electromotive force across a resistance R 



u:. 



di 
is Ri and that across an inductance L is L -tt' The e.m.f. impressed 

upon a circuit containing a resistance R and an inductance L in series 

is 

e B ^ sin (cd), 

E and w being constants. Find the current at time t U i — when 
t = 0. 
. 4. Two circuits of resistance JBi, Rt and inductance Li, Lt respec- 
y tively are connected in parallel between the mains of a transmission 

Hne. If the total current they receive is 

• 

i B / sin («0, 

find the ciurent in each circuit and the e.m.f. between the mains as- 
suming that both currents are zero when ^ = 0. 

5. A cylinder containing gas is rotated with constant angular velocity 
a> about its axis. Ultimately the mass of gas will rotate like a rigid 
body. Assuming Boyle's law and taking account of centrifugal force, 
find the law connecting the pressure of the gas and the distance from 
the axis. Would a similar law be obtained, if instead of a cylinder, 
the container had any other shape? 

6. In a chemical reaction a substance c decomposes into two sub- 
stances X and y, the rate at which these products are formed being pro- 
portional to the amount of- c present. If at the beginning c = 1, 
X — 0, y == Of and at the end of 1 hr. c = i, a; = J, y = j, find x and 
y as functions of the time. 

7. In a certain chemical reaction, 1 mol. of 2^ is produced for each 

mol. of X consumed, the rate being proportional to the amount of x 

present, and at the same time y by & reverse reaction is converted into 

X at a rate proportional to the amount of y present. Chemical analysis 

showed 

< = 0, 3, 00 

x = 10, 6, 5.5 

y = 0, 4, 4.6. 

Find X and y as functions of t 
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8. Radioactive substances are transformed at a rate proportional 
to the amount of substance present. Through the decomposition of 
a mol. of Ra B a mol of i2a C is produced, the rate being such that i 
the Ra B disappears in 27 minutes. Similarly the Ra C decomposes 
at a rate such that i of it is lost in 19.5 min. If initially 1 mol. of 
RaBia present, find the amount of Ra B and RaC At the end of 1 hr. 

9. When light passes from a medium of refractive index m to one of 
index ti\ 

n' sin0 
M sm0' 

and $' being the angles which the incident and refracted ray make 
with the normal to the surface of separation. According to Einstein's 
theory, the gravitational field of the sim deflects a ray of light as if 
it had a refractive index. 

r 

where a is constant and r the distance from the center of the sun. 
Find the path described by the ray. 

10. Suppose bacteria grow at a rate proportional to the nimiber 
present but that they produce toxines which destroy them at a rate 
proportional to the nimiber of bacteria and to the amount of toxin. 
Suppose further that the rate of production of toxin is proportional to 
the number of bacteria. Show that the number increases to a maxi- 
mum and then decreases to zero and at time t is given by 

4Af 



N = 



(e« + e-*')2 



where M is the maximum number and t is measured from the time 
when the number is a maximum. 



CHAPTER III 
SPECIAL TYPES OF SECOND ORDER EQUATIONS 

18. Equations Immediately Integrable. — An equation 
of the form 

g = /(x) (18) 

can be solved directly by two integrations. The first 
integration gives 

A second integration gives the general solution in the form 

y = / ( I f{xydxjdx+ cjx +-C2. 

19. Equations not Containing y, — An equation not con- 
taining y can be solved for the second derivative and so 
reduced to the form 



g = ff. ^\ 



Take as new variable 



p-t' (i«b) 



Then 

(Py _dp 
dx^ dx 

and so (19a) can be written 

33 
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This is a first order equation whose solution has the form 

p ^ F(x, ci), . 

where Ci is the constant of integration. 
Substituting ~ for p, we have 

whence 



y =JF(x,Ci)dx + C2. 



Exampk. (l + x)g + ^=0. 

dv 
Substituting p for -^ we have 

This is an exact equation with solution 

(1 + x)p = Ci. 

di/ 
Replacing p hy -^ and separating the variables, we have 

T Ci dx ^ 

whence 

2/ = Ci In (1 + x) + C2. 

20. Equations not Containing x. — An equation not con- 
taining x can be reduced to the form 

§■=/(!/. P)' (20a) 

Substitute 

| = P (20b) 
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and write the second derivative in the form 

da? dx dydx dy* 
These substitutions bring (20a) to the form 

P^ = f(y.P). 

This is a first order equation which can be solved for p. 

dv 
Replacing P^Y jf- the result is a first order equation which 

can be solved by a separation of the variables. 

Substituting 

dx~^' d3? '^dy 



the equation becomes 



Dividing by p, 



yv-£='y^ + T^' 



yfy = y' + v. 



The solution of this equation is 



whence 



-y-y + c., 



| = .(y + cO 



and 



-/ 



_^=lln^ + C. (20d) 

y{y + ci) «i y + ci 
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In solving this problem 'we divided both sides of the 
equation by p. This is allowable since 

p = 
gives 

y = c (20e) 

a solution which contains only one constant of integration 
whereas we are determining the general solution with two 
constants of integration. It is to be noted that (20d) is a 
solution of the original differential equation and that it 
cannot be obtained by giving special values to the con- 
stants in (20c). Such a solution is called singvlar and in 
some problems might be important. 

EXERCISES 

Solve the following differential equations: 

5. gt*^=0- 



6. 






9. (. + l)g-(. + »| + . + 2-0. 
>»■ »g - ' + (I)' ■ 



12. 
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21. Deflection of Beams. — When a beam is bent by 
vertical forces as shown in Fig. 21a, the fibers in the upper 
part are stretched and those in the lower part compressed. 
There is then a neutral curve. AB along which they are 
neither stretched nor compressed. It is seen from the 
figure that the amount ds a 
fiber of natural length s is 
stretched or compressed is 
given by the equation 

ds __ s 

where z is its distance from 

the neutral curve and 22^ is 

the radius of the circle in 

which the neutral curve is bent. By Hooke's law the tension 

in the fibers between z and z + dz is 

T = Ewdz^^ = Ewdz.^f 

s R 

where E is the stretch modulus of elasticity of the ma- 
terial and w the width of the beam. Since there is no re- 
sultant force along the beam 

/ Ew dz>r^==-^jwzdz = 

showing that the neutral curve passes through the center 
of gravity of each cross section. The moment of the total 
stress about an axis PQ perpendicular to the beam is 



M 



RJ 



wz^ dz = 



R 



(21a) 



when / is the moment of inertia of the section area about 
that axis. 
Let (x, y) be the coordinates of P, x being measured 
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along the beam. Since the curvature is usually very small, 
the slope -p is nearly zero and so 

1 dx^ ^^ 



['+(!)■] 




X 




Fig. 21b. 



approximately. Hence (21a) can be written 

S7g=ilf. (21b) 

M, called the bending moment, is the moment about P 
of all forces on one side of P, those acting upward producing 
positive moments, those downward negative. To find the 
Y curve in which a beam bends 

we determine ilf as a func- 

I [ } I I I { tion of X and solve the 

■ ^ "Ll ^ differential equation (21b). 

In the problems considered 
here the beam is of uniform 
material and has a constant 
cross section. Hence E and / are constants. 

Example. Find the deflection of a beam of length 2 Z, 
supported at its ends and loaded with a weight w per unit 
length. 

Take the origin at the center of the beam. The total 
weight supported being 2 Iw, the upward thrust of the sup- 
port at each end is Iw, Consider forces on the right of the 
point P. At the end is an upward thrust Iw. Its moment 

about P is 

wl (I — x). 

The only other force on the section PA is the load w (I -- x) 
between P and A. Since the load is uniformly distributed, 
this can be considered as acting through its center of gravity 
at distance i (I — x) from P. Its moment is then 

— w (I — x) • ^ (I - x), 



Chap. Ill SECOND ORDER EQUATIONS 39 

the negative sign being used because the force is downward. 
The total moment about P is 

M ^ wiq, - x) -"^(J, - xf = \w Q* - a?). 
The differential equation is therefore 

A first integration gives 

Since the beam is horizontal at the center 

andsoci = 0« A second integration gives 

Since we have taken the origin on the curve ct = 0. The 
equation of the elastic curve is therefore 






The maximum departure of the beam from a straight Une is 

at the end, given by placing x = Z. This is ordinarily- 
called the deflection of the beam. 

22. EquiKbrium of a Cable. — Let ABC (Fig. 22) be a 
perfectly flexible cable fastened at A and C and loaded in 
any way. Take the x-axis horizontal and the 2/-axis through 
the lowest point B on the curve. Consider the part of the 
cable between B and the variable point P (x, y) on the 
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curve. This is in equilibrium under the action of three 
forces: 

(1) A horizontal tension H at B exerted by the section 
AB of the cable. 

(2) A tension T along the tangent at P due to the part 

PC of the cable. 

(3) A downward force 
equal to the load W on 
the part BP of the cable. 

The total component of 
force toward the right must 
equal that toward the left, 
and the component of force 
acting upward must equal 
that acting downward. 




Fig. 22. 



Hence 

r cos = £r, 

Tsin* = W. 
Dividing the second equation by the first, we get 

x' dy W 

tan = -;- = 77 • 
ax H 



(22a) 



(22b) 



Since H is constant, if TT is a known function of x, this can 
be integrated at once. In some cases W is not known but 
its derivative can be easily determined. In that case, 
differentiation gives 

(Py 1 dW ,^ V 

from which we can determine the equation of the curve. 

The answer contains the constant H, This can be de- 
termined by substituting the known coordinates of B and 
the ends of the cable. 

23. Motion of a Particle in a Straight Line. — If /^ is the 
resultant of all forces acting on a particle of mass m, its 
acceleration a is given by the equation 

F ^ma. . (23a) 
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If the particle moves along a straight line and 8 is its 
distance from a fixed point of the line, its velocity is 

v^%, (23b) 



and its acceleration is 



dv d^8 /oo \ 



In using these formulas we can measure 8 in either direc 
tion along the line but must then consider as positive the 
direction in which 8 increases. The quantities F, v, and a 
are positive or negative according as they point in the 
positive or negative direction thus defined. 

Example. When a body sinks slowly in a Uquid the 
resistance is approximately proportional to the velocity. 
If the particle starts from rest, find its motion. 

Let s, considered positive downward, be the distance 
the body sinks in t seconds. If m is its mass and g the ac- 
celeration of gravity, the force of gravity is 

mg, 

which is positive since the force is downward. The re- 
sistance acting upward and being proportional to the ve- 
locity is 

— h). 

The total force acting on the body is then 

F = mg — kv. 

Hence equation (23a) is 

, dv 

mg — kv =^ ma — m-r:' 

at 

Separating the variables and integrating, we get 

k 
In (mg — fct;)= t + c 
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Since t; = wheu i == 

In {mg) = c. 

Subtracting from the preceding equation and solving for 
mg — ku, 

t 

mg " kv == mg e~m\ 

mgt — ks = j^ e m* + c. 

K 



Since s = when i = 0, 






Substituting this value and solving for s, 

24* Motion of the Center of Gravity, — If M is the total 
mass of a body or system of bodies and F the resultant of all 
the forces appUed to it, the equation 

F = Ma (24) 

determines the acceleration a of the center of gravity. 
That is, the center of gravity moves as if the whole mass 
were concentrated at that point and all the forces applied 
there. 

If all parts of the body move in the same direction with 
the same velocity, this equation determines the acceleration 
of any point of the body. If the parts of a complex system, 
such as a chain, all move along the same path with the 
same speed, and F is the component of force along the path, 
(24) gives the acceleration along the path. 
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26. Motion in a Plane. — When a particle of mass m 
moves in a plane or in space^ its acceleration still satisfies 
the equation 

F = ma 

but the quantities F and a are vectors, that is, have direc- 
tion as well as magnitude. To obtain an equation whose 
terms are nmnbers we project on any line. If the com- 
ponent of F along any line is F^ and the component of a 
along the same line is a^, then 

Fx = TTULx* 

Suppose the particle moves in a plane. Let (x, y) be its 
rectangular coordinates. The components of acceleration 
along the axes are 



a* = 






Oy = 



_d^y 



de 



If Fx and Fy are the components of force along the axes, 
the motion of the particle can then be determined by solving 
the differential equations 



cPx „ 

In problems where the force 
acts along the line joining the 
particle to a fixed point it may 
be more convenient to use 
polar coordinates with that 
point as origin. When the 
particle is at P its accelera- 
tion is resolved into a com- 
ponent Or along OP and a 
component Oe perpendicular to 
OP. These are 






(25a) 




dV /dev 



a 



B 



Fig. 25a. 



rdt\ dtj 



(25b) 
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If the force acting on the particle has the component 
Ft along OP and F^ perpendicular to OPy the motion of the 
particle can be found by solving the equations 



mor = Ftj 

ma^ = Fg. 
To prove equations (25b) write 



(25c) 



X ^ rco8d,y == rsind 



and calculate 









in terms of r and 6, These are accelerations along OX 
and OY, To obtain Or project a^ and Oy on OP and add the 
results. Similarly, to obtain o^ project on the line per- 
pendicular to OP. 
Example 1. When an electron of charge e moves with 

velocity v in a field of mag- 
netic intensity H, it is pushed 
sidewise with a force 



e 




F = - rff sin 
c 



Fig. 25b. 



where c is the velocity of light 
and the angle between v 
and H. The force is perpen- 
dicular to both V and H and 
so directed that a right hand rotation about F through the 
angle 4> carries v into H. 

If H is constant and the electron starts with velocity 
Vo in a plane perpendicular to H, find the path de- 
scribed. 
Let the electron start at the origin (Fig. 25c) with initial 



Chap. Ill 



SECOND ORDER EQUATIONS 



45 



velocity along the x-axis. Since there is no component of 
force perpendicular to the xy-plane, it will continue to move 
in that plane. Hence <t> = 90° 
and 

c 
The components of t; are 



dt 



dy 
dt 



Since F is perpendicular to 
V and numerically equal to 

- vHf its components are 

e„dy 
c dt 




Fig. 25c. 



e Yjdx 
c dt 



the algebraic signs being determined by inspection of the 
figure. The equations of motion are therefore 



m 



d^x __ _ eHdy 
dp " ~7di 

d^y _ eHdx 
c dt 



"*d?= 



These equations can be immediately integrated giving 



dx 



m-TT = 



dt 
dy 



eH , 





eH 



m-ir = 



dt 



Xs 



the constai^ts of integration being determined so that 

^ r, dx dy g^ 
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when < ~ 0. Dividing the second equation by the first 



dy eHx 



dx cmvo — eHy 

Clearing fractions and integrating, we get 

eH (a^ + y^) — 2 cmvo y = 0, 

the constant being zero since y == when a? = 0. The 

electron therefore describes a 
circle of radius 

cnwo 

Example 2. A particle of 
mass m is attracted toward 
the origin with the force 




If it starts from the point (a, 0) with velocity vo> - per- 
pendicular to the X-axis, find the motion. 
Using polar coordinates, the differential equations are 



m 



Idt* \dtj J 
md/,de\ _ Q 
r dt\ dt) "• 






The second equation can be inunediately integrated giving 



.de 



At the start r = a and 



dB 
^dt 



= f^b. 
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Hence Ci = avo and 

d$ _^ 

Substituting this value in the first equation 



dp r« 



dv 
multiplying by ^rr dt, 



dt 



whence 



dr dV j^ aW — fc^ j 

5J?* p— * 



(D'-(.v-».,(i-i) 



dr 



the constant of integration being determined so that jt = 

at 

when r ^ a. Dividing by 



\dt) r* 



we obtain 



(dry ^ ^ 
W 



aW 



r^ (r^ — a^). 



Solving this equation and determining the constant so that 

r =^ a when ^ = 0, we finally obtain the equation of the 

curve in the form 

rVaW - k^l 
r = a sec 6 • 

L CiVo J 

26. Rotation about a Fixed Axis. — Let the body (Fig. 
26) rotate about the axis through perpendicular to the 
plane AOB. The position of the body can be determined 
by the angle 6 between.the fixed line OA and the Une OB 
fixed in the body. The angular velocity of rotation is 
then 

« = I (26a) 
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and its f^^n giilar acceleration 



d<a cPB 



a = TT = 



dt d^ 



(26b) 



Let F be a force in the plane OAB applied to the body at 
P. The torque about the axis of rotation due to this force 

is 

r = FZ, 

where I is the p«pendicular distance from O to the line FP. 
The torque is positive when the force F tends to increase 0, 

The torque and angular 
^ acceleration satisfy the equa- 
tion 
B r = /a, (26c) 




where / is the moment of 
A inertia of the body about the 
axis of rotation and T is the 
YiQ, 26. total torque about that axis of 

all forces acting on the body. 
This is analogous to (24), the moment of inertia correspond- 
ing to mass and torque corresponding to force. 

27. Combined Translation and Rotation. — Consider the 
motion of a rigid body of mass M whose center of gravity 
moves in a fixed plane and which at the same time rotates 
about an axis perpendicular to that plane. In the problems 
solved here the axis is always an axis of symmetry and the 
forces acting on the body lie in the fixed plane. The axis 
of rotation will then remain perpendicular to that plane. 
The motion of the center of gravity is determined by the 
vector equation 



F = Ma 
and the rotation by the equation 



(27a) 



(27b) 
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In case of a system consisting of two or more rigid 
bodies moving as just described there is an equation of 
the form (27a) and one of the form 
(27b) for each body. The force 
F is in each case the resultant of 
all forces acting on the body. It 
may not be possible to determine 
at once the force that two bodies 
of the system exert on each other. 
Such a force can be represented 
by a letter. It will be found that 
there are enough equations to 
determine these imknown forces as 
well as the accelerations. 
. Example 1. A cyUnder of mass 
M and radius r rotates about its 
axis. A cord wrapped aroimd the fig. 27a. 

cylinder is attached to a mass m 

which drops vertically. If the cyUnder starts from rest, 
find the angle turned through in t seconds. 

Let F be the tension in the cord. The torque about the 
axis of the cylinder is 

T == Fr. 

The moment of inertia of a cyUnder is 

I = Imt^. 

Equation (26c) is then ' 

. ^Mt^ a = Fr 

or 

^Mra = F. (1) 

The forces acting on m are F acting upward and the 
force of gravity mg acting downward. Its -acceleration 
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then satisfies the equation 

ma =ing — F. (2) 

If is the angle the cylinder turns through in t seconds 
and 8 the distance m falls 

s = rO 

and so 

(Ps cPB 



or 



dt^ ^dfi 



a = ra 



(3) 



By solving (1), (2), (3) simultaneously, we obtain 

2 mg d^e 



a 



Hence 



e = 



(Af + 2 m) r d? 
(Af + 2m)r* 



Example 2. A sphere of mass M and radius r rolls down 
a plane which makes an angle with the horizontal. If 

the coefficient of friction is n 
and the inclination is so great 
that the sphere slides, find its 
angular acceleration about the 
horizontal axis through its 
center of gravity. 

The force of gravity on the 
sphere can be resolved into 
Fig. 27b. ^^^ components 

F = Mg sin 0, P = Mg cos 

parallel and perpendicular to the plane. The force of fric- 
tion is 

fiP =^ n Mg cos 4>. 

The torque due to this force is 

r = M Mg r cos <f>. 




Chap. Ill SECOND ORDER EQUATIO^JS 61 

The moment of inertia of a sphere is ; 

I = %Mr^. 
5 

Hence 

2 

/i Af gr cos <^ == -zMr^ a 



and so 



5 un cos 
2 r 



PROBLEMS 



1. A beam of length 2 Z is supported at its ends and loaded with a 
weight W at the middle. Find the deflection. 

2. Find the deflection of a cantilever beam of length Z, held horizon- 
tal at one end, and loaded with a weight W at the other end. 

3. Find the deflection of a cantilever beam fixed at one end and 
loaded with a weight w per imit length. 

4. Find the deflection of a beam supported at both ends and at the 
middle point, loaded with a weight w per unit length. 

5. Find the deflection of a cantilever beam fixed at one end, sup- 
ported at the other and loaded with a weight w per imit length. 

6. Find the deflection of a beam fixed at both ends and loaded with 
a weight w per unit length. 

7. ConJsider a vertical column fixed at the base, of length Z, and sup- 
porting a weight P. Suppose the weight causes the upper end to be 
displaced the amount a from the vertical. Calculate the bending 
moment and determine the curve in which the column bends. By 
substituting the coordinates of the upper end show that the maximum 
load the column can support is 



- («)■ 



EL 



8. Suppose the ends of the column are roimded but are held in the 
same vertical line. Find the maximum load the coliunn can support. 

9. The cable of a suspension bridge supports a bridge of weight w 
per unit horizontal distance. Neglecting the weight of the cable find 
the curve in which it hangs. 

10. A series of rods of varying length but the same diameter are 
hung along a cord. The horizontal distances between consecutive 
rods are equal and their bottoms are in a straight line. Assuming 
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that they are so close together that the load can be considered con- 
tinuous, find the curve formed by the cord. 

11. A cable is supported by its ends and hangs under its own weight. 
Find the curve in which it hangs. 

12. A telegraph wire weighs 173 lbs. per mile. If the poles are 400 
ft. apart and the wire sags 10 ft. at the middle, find the tension at the 
lowest point of the wire. 

13. An arch of a masonry bridge supports a horizontal roadbed and 
is so constructed that the resultant stress at each point of the arch 
due to the material above is a compression along the tangent. Find 
the shape of the arch. 

14. A particle of mass m moves in a straight line toward a center of 
force which attracts with the magnitude 



where r is the distance from the center. If the particle starts from 
rest at the distance a, find the time required to reach the center. 

16. A motor boat weighing 1000 lbs. is moving in a straight line with 
a velocity of 50 ft./sec. when the motor is shut off. If the resistance 
of the water is proportional to the velocity of the boat and is equal to 
10 lbs. when the velocity is 1 ft./sec., how far wiU the boat move before 
the velocity is reduced to 25 ft./sec. How long will be required for 
this reduction in velocity to take place? 

16. A particle of mass m moves toward a fixed center of force which 
repels with a force khn times its distance from the center. If it starts 
from the distance a with velocity koj show that it will continually 
approach but never reach the center. 

17. Find the velocity acquired by a body falling from an indefinitely 
great distance to the earth, assuming that the force of attraction varies 
inversely as the square of the distance from the center of the earth. 

18. Find the time required for a body to fall to the earth from a 
distance equal to that of the moon. Take the radius of the earth as 
4000 mi. and the distance from the center of the earth to the moon as 
240,000 miles. 

19. If a hole were bored through the center of the earth, a body 
, falling in it would be attracted toward the center with a force propor- 
tional to the distance from the center. Find the time required to 
fall through. 

20. A body slides down a rough inclined plane. If the inclination 
of the plane is a and the coefficient of friction is fi, determine the motion 
if the particle starts from rest. 
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21. Assume the resistance of the air proportional to the square of • 
the velocity. If the velocity of a falling body is observed to approach 
the limiting value 216 ft./sec. and the body starts from rest, find the 
motion. 

22. A particle is projected vertically upward with velocity vq. As- 
suming that the resistance of the air is A; times the square of the velocity, 
find the velocity with which it returns to the earth. 

23. A chain 6 ft. long starts with 1 ft. of its length hanging over the 
edge of a smooth table. Neglecting friction, find the time required to 
slide off. 

24. A chain hangs over a smooth peg, 8 ft. of its length being on one 
side and 10 ft. on the other. Find the time required to sUde off. 

26. Solve the preceding problem if the force of friction is equal to 
the weight of 1 ft. of the chain. 

^6. A projectile is fired with a velocity of 2500 ft./sec. in a direction 
making 45° with the horizontal. Find the highest point reached and 
the point where it strikes the ground. 

27. A projectile is fired with velocity vo at an angle of elevation a. 
If the resistance of the air is kv, where v is the velocity and A; is con- 
stant, find the equations of motion. 

28. A particle is attracted toward the origin with a force propor- 
tional to the distance. If it starts from the point (a, o) with velocity 
t^ perpendicular to the x-axis, find the path described. 

29. Solve the preceding problem if the particle is repelled with a 
force proportional to the distance. 

30. An electron moves in a magnetic field of intensity H. If it M 
starts with velocity vo in a direction making the angle a with H, find ^ 
the path described (See Example 1, page 44). y^ 

31. Prove equations (25b). ^' 

32. Determine the orbit of a planet of mass m assimiing that it is 
attracted toward the sun with the force 

km 

-^' 

where r is the distance from the sun. Let ro be its distance and vo 
its velocity when nearest the sim. 

S3. Find the orbit of a comet. Let ro be its least distance from the 
Sim and assume that its velocity at an infinite distance is zero. 

34. A particle of mass m is attracted toward the origin with the 
force 

^.^.^ • 

2* 
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If it starts from the point {a, 0) with velocity 

k 



vo = 



a*y/2 



perpendicular to the x-^axis, find the path described. 

36. A circular disk of radius a submerged in oil rotates about the 
perpendicular axis through its center. Assume the frictional resistance 
per imit area at each point of the disk to be A»;, where v is the velocity 
of that point and k is constant. If the disk is started with angular 
velocity c^and the torque due to friction at the bearings is a constant 
K, find the motion. 

36. A ball of radius r rolls without slipping down a plane. If 4> is 
the angle of inclination of the plane and the ball starts from rest, find 
the distance its center moves in t seconds. 

37. A billiard ball is started with velocity vo not rotating. If the 
coefficient of friction between the ball and table is n, find the motion. 

88. A cylinder of mass M and radius r rolls on the top of a table. 
A cord wrapped around the cylinder passes horizontally over a fixed 
pulley and is attached to a weight m which drops vertically. Find the 
motion of the cylinder. 

39. Solve the preceding problem if the cord is attached to the axis 
of the cylinder. 

40. A wedge shaped block of mass M and 45° angle slides on a 
smooth table. A mass m sUdes on its surface. If both start from 
rest, find the motion. 



CHAPTER IV 
/ 

LINEAR EQUATIONS WITH CONSTANT 

COEFFICIENTS 

28. Equations of the n-th Order. — The solution of an 
equation of the n-th order in two variables involves n in- 
tegrations. The general solution therefore contains n con- 
stants of integration. 

Constants of integration are called independent if they 
occur in the solution in such a way that it is not possible 
to replace a function of two or more of them by a single 
constant and so reduce the number of constants. Thus 

2/ = ciar^ + C2 + C8 (28) 

appears to contain 3 constants. We can however take 

C2 + Ci == c^ 

and so obtain 

y = Cix^ + c. 

It is not possible to further reduce the number. Hence 
(28) contains two independent constants. 
In case of a differential equation of the form 






where / is an algebraic function, it can be shown that there 
is only one solution containing n independent constants of 
integration. It is called the general solution. 

A solution obtained by giving particular values to the 
constants in the general solution is called a particular 
solution. Some differential equations have solutions con- 
taining less than n constants of integration which are not 

55 
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particular solutions. Such solutions are called singular. 
They are mainly of mathematical interest and so will not 
be further considered here. 

29. Linear Equations with Constant Coefficients. — A dif- 
ferential equation of the form 

is called a linear equation. By this it is meant that the 
equation is of first degree in one of the variables (y in this 
case) and its derivatives. If the coefficients ai, 02, . . , 
On are constants, it is called a Unear equation with constant 
coefficients. For practical appUcations this is one of the 
most important types. 

In discussing these equations we shall find it convenient 

to represent the operation -y- by D. Then 

az 

Equation (29a) can be written 

' (D» + aiD«-^+ . . '+a^,D + an)y=f{x). (29b) 

This signifies that if the operation 

D» + aiZ)»-^ + • • • + a^.Z) + a« (29c) 

is performed on y, the result is /(x). The operation con- 
sists in differentiating y, n times, n — 1 times, etc., mul- 
tipljdng the results by 1, ai, 02, etc., and adding. 
With the differential equation is associated an algebraic 

equation 

r» + aiT^" + . . . + an-,r + a„ = (29d) 

having the same coefficients ai, 02, etc. as (29a) but with 
right-hand member zero. If the roots of this atixiliary 
equation are n, r2, • • • , fn, the polynomial (29c) can be 
written 

, . (D-n) iP'-r^)^ ^ *(D-rn) 
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and so (29a) has the form 

(D - n) (D - rj) . . . (D - rn) y = /(x). (29e) 
If we operate on y with D — ri, we get 

iD-n)y = ^£- ny. 

/ 

Operating on this with D — r2 we get 

(D - n) (£> - rO y = (D - r,) (g - r,j/) 

= ^ - (n + r2) 1/ + rir2y. 

The same result is obtained if we operate on y with 

{D - n) (D - r2) = Z)2 - (n + rj) D + rir2. 

Similarly, if we operate in succession with *the factors 
(D — n), (D — r^, etc., in any order whatever we get the 
same result that we should get by operating directly with the 
product (29c). It should be noted that this is true only when 
ri, r2, etc., are constant. If the r's are variable it is not in 
general true that 

(D - n) (D - rj) 2/ = (D - r2) (D - n) y. 

30. Equation with Right-hand Member Zero. — To solve 
the equation 

(D» + aiD»-' + • • • + a,^,D + On) 2/ = 0, (30a) 

factor the symbolic operator and so reduce it to the form 

(D - n) (D - r2) • • • (D - r„) y = 0. 

The value 

y = Ci6^ 

is a solution. For 

(D — ri) cic^' = Cine"* — ncie"* = 0, 
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and the equation can be written 

{D-n) . . . {D-rn) • (D-n) y=(D-r2) • • • (Z>-r„) • 0=0. 

Similarly 

y = c^, y = ct^f etc., 

are solutions. Finally 

y = Ci6"* + 02^'"* + • • • + CneV (30b) 

is a solution; for the result of operating on y is the sum of 
the results of operating on Cic'"*, c^^*', etc., each of which 
is zero. 

If the roots n, rj, • • • , r» are all different, (30b) contains 
n independent constants and so is the general solution of 
(30a). 

If, however, two roots n, r2 are equal 

CieT^ + C2e''*' = (ci + Ci)e''»' 

contains only one constant Ci + C2 and (30b) contains less 
than n independent constants. In this casej however, 
Qc^^ is a solution; for 

(D — n) xe*"^* = rixe*"'* + c'^* — nxe'^* == c'^* 

and so 

(D - ri) (D - Ti) xe'»* = (D - ri)^ xe^^^ = (D - n) e'»* = 0. 

If then 

r2 = n 

the part of the solution corresponding to these two roots is 

(ci + C2x)e'**. 

More generally, if 

fi = r2 = • • • = rm, 

the part of the solution corresponding to these m roots is 

{ci + dpc + csx^ + ' ' ' + Cnai^-')e'^ (30c) 
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If the coefficients ai,ai,* • • On are real, imaginary roots 
occur in pairs 

ri = a + /S V"^, r2 = a - /S V^. 
The terms 

are then imaginary but in most problems their sum is real. 
They can be replaced by two other terms that do not have 
this imaginary appearance. Using the values of n and r2 
we have 

(D - n) (D - r2) = (D - a)« + 0^. 

By performing the differentiations it can easily be shown 
that 

[(D - ay + /?] . e«*sin/Sx = 0, 

[(D - ay + /32] . e«* cos Px = 0. 

Therefore 

^ [ci cos j3x + Ci sin /Sx] (30d) 

is a solution. This fimction in which a and j9 are real can 
therefore be used as the part of the solution corresponding 
to two imaginary roots 

To solve the differential equation 

(D» + aiD»-^ + • • • + a„_,D + On) 2/ = 
let n, r2, • • • , r« be the roots of the avxiliary equation 

r" + air«"> + •••*+ an-^r + On = 0. 
7/ iA6S6 roots axe real and different the solviion is 

y = cie^ + c^"^ + • • • + Cne^n^. 

// m of the roots ri, r2, • • • r„ are equals the correspondr 
ing part of the sohdion is 

(ci + Cax + CiX* + . . . + c„^«i-i) ^. 
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The part of the solution corresponding to two imaginary 
roots r =^ azL P V— 1 is 

^ [ci coe jSx + c» sin jSx]. 

This is equivalent to 

(Z)« - Z> - 2) y = 0. 
The roots of the auxiliary equation 

r* - r - 2 = 
are — 1 and 2. Hence the solution is 

y = Ci6"* + o^ 

ExamvU2. g + g-5|+3y = 0. 

The roots of the auxiliary equation 

are 1, 1, — 3. The part of the solution corresponding to 
the two roots equal to 1 is 

(ci + Cax)c*. 
Hence 

y = (Ci + C^) €F + C8€"^. 

Example 3. {D^ + 2D + 2) y = 0. 
The roots of the auxiliary equation are 

- 1 zb V^T. 

Therefore 

a = - 1, /3 = 1 in (30d) 

and 

y = e-« [ci cos x + Cs sin x]. 

31. Equation with Right-hand Member a Function of x. — 

Let 2/ = w be the general solnJtion of the equation 

(D« + aiZ)»-^ + • • • + a«_,D + On) 2/ = 
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and lety — vhe any solution of the equation 

<D» + aiD»-« + - • * + a^J) + an)y^f{x) (31) 
then 

is a solution of (31) ; for when the operation 

2)»+ aiD»-^ + • • • + a^,Z> + a, 

is performed on u it gives zero and when it is performed 
on V it gives / (x). Furthermore u + v contains n ar- 
bitrary constants. Hence it is the general solution of (31). 

The part u is called the complimentary function^ v the 
particular integral. To solve an equation of the form (31) 
we first solve the equation with right-hand member zero 
and then add to the result any solution of (31). 

A particular integral can often be found by inspection. 
If not the geiieral form of the integral can be determined 
by the following rules: 

1. Uf{x) = ox" + fex"-' +••• + ?>, 

assume 

y = i4x» + Bx"-* + . • • + P, 

but, if occurs m times as a root of the auxiliary equation, 
assume ' ' 

y = aJ« [i4x« + Bx"-*' + • • • + P]. 

2. If/(x) = ce^, assume 

y = Aef^, 

but, if a occurs m times as a root in the auxiUary equation, 
assume 

y = ilx" €f". 

3. Uf{x) = a cos iSx + 6 sin fix, assume 

y = A cos fix + Bmi fix, 
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but, if COS fix and sin fix occur in the complementary func- 
tion, assmne 

y = x[Acoafix + Bmi fix]. 

4. If / (x) = a^ cos fix + b^ sin fix, assume 

y = AeT" cos fix + BeT" sin fix, 

but, if ef* cos fix and ^ sin j9x occur in the complementary 
function, assume 

y = xef^ [il cos jSx + J? sin fix]. 

If / (x) contains terms of different types, take for y the 
sum of the corresponding expressions. Substitute the as- 
sumed value of y in the differential equation and determine 
the constants A, B, C, etc., so that the equation is satisfied. 

The general principle in the above rules is to express y 
as a Unear function of all the distinct kinds of functions in 
/ (x) and its derivatives of all orders. The exceptions to 
the various rules occur when some of the terms in the as- 
sumed value of y occiu* in the complementary function. 

Example 1. j^ + 42/ = 2x + 3. 



dxi 



A particular integral is evidently 



The solution of 



j/ = i(2x + 3). 



g+4„-0 



IS 



2/ = Ci COS 2 X + Ci sin 2 X. 
Hence the solution of the original equation is 

y = Ci cos 2 X + C2 sin 2 X + 2 (2 X + 3). 

ExampU 2. (D« + 3 D + 2) y = 2 +<f 
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Assume 

y == A + B^. 

Substituting this value for y 

2 A + 6Be* = 2 + e*. 
Hence 



and 



Hence 



2il = 2, 6B = 1 



A + B^=l + ^^. 



y = 1 + ge* + cie-*+ c^""^. 



Examples. ^ + g=:r?. 

The roots of the auxiUary equation are 0, 0, — 1. 
Since is twice a root, we assume 

y ^ x^ {Aj? + Bx + c). 
Substituting this value, 

l2A7?+{2AA^-6B)x + &B + 2c = 7?. 

Consequently 

12A = I, 24A + 6B = 0, 6B + 2c = 0, 

whence 

The solution is 

y = j2^ "■ 3^ + ^ + Ci + C2X + cse-*. 

3J2. Simultaneous Linear Equations. — We consider only 
linear equations with constant coefficients, containing one 
independent variable and as many dependent variables as 
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equations. All but one of the dependent variables can be 
eliminated by a process analogous to that used in solving 
linear algebraic equations. The one remaining dependent 
variable is the solution of a Unear equation. Its value can 
be found and the other dependent variables can then be 
determined by substituting in the preceding equations. 

If possible the work should be so arranged that after the 
first variable is found the others can be determined without 
integration. If integration is used in determining these 
later variables, the constants of integration may not all be 
arbitrary. It is then necessary to substitute the values 
found in the differential equations to determine the relations 
between the constants. 

Example. -^ — Sx — y —. ef 

at 

Using D for t;» these equations can be written 

(Z>2 - 3) a; - 2/ = e* 
'-2X + Dy = 

Multiplying the first equation by 2 and the second by D^ — 3, 
we have 

2(Z>2 -3)x-22/ = 2^, 
- 2 (Z)2 _ 3) X + (D2 ^s)Dy = 0. 

Adding, we get 

(Z)« - 3 D - 2) y = 2 e*. 

This equation, containing only one dependent variable, 
can be solved for y giving 

Substituting this value of y in the second equation, we find 
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EXERCISES 

db^ dx ' - 

>>4. g + y-O. 

5 ^-2^-3^ = 
-^ ■ das* das* dx 

10. g + ^-2-x. 

IS. ^ — y •■ sin «. 

18 ^_*?«a! 
da^ dx 

15. g-«V-e«. 

16. g + o?y-ooBOK. 

17. g^y-a^-L 

19. ^ - 9 y = ««« cos «. 
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dx i t dy , . 

dx n dy 

*^- S^ + 5i + *"^' di + d^'^- 

33. Vibrating Systems. — When a body is given a slight 
displacement from a position of stable equiUbrium and re- 
leased, it vibrates about the position of equihbrium. 

If the resistance is neglected the differential equation of 
motion usually has the form 

The solution of this equation can be written 

X = il cos (fc< + e), (33a) 

when A and are constants. The motion is called har- 
monic. A is called the amplitude and 6 the phase angle. 
In a complete vibration (across and back) the angle 

kt + e 
increases 360®, or 27r. The time of vibration is therefore 

If the resistance is proportional to the velocity, the dif- 
ferential equation has the form 

cPx . dx , , ^ 

The solution of this equation has the form 

X = Ae-^' cos (m + 6). (33b) 
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If at time t, the body is at the end of a swing, 

is the amplitude of that swing. The amplitude thus de- 
creases with the time. 

In some cases the differential equation has the form 



dt^ 



+ fc2 sin ^ = 0. 



This is not a linear equation. It becomes linear, however, 
if we replace sin d by 6. For small angles this is a good 
approximation. Thus, when the 
angle is 5° 

sin ^ =^ .08716, 
B = .08727, 

showing that the error is less than 
1 part in 800. Even for angles of 
10° the error is only about \ per 
cent. 

Exampk. A pendulum, consist- 
ing of a particle of mass m sup- 
ported by a string of length I, 
swings in a medimn which resists 
with a force proportional to the velocity. Find the time 
of vibration. 

The torque about the point of suspension due to the 

weight of m is 

— mg sin ^ 

the negative sign being used because when d is positive the 
torque tends to decrease 6. The velocity of m is 




Fig. 33. 



,de 



The resistance is then 



kl 



dd 
dt 



I 
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A 

The torque diig^to this resistance is 

1.7 ^ I 

the negative sign being used since when -^ is positive the 

torque is negative. The moment of inertia of m about 

is mP. 

The equation 

J— = r 

is therefore 

at at ' 

Replacing sin 6 by 6, this becomes 

The solution of this equation has the form 



^ = il e""2m cos 



The period of vibration is therefore 

27r 



r = 



V Z 4m^ 



2 



PROBLEMS 

1. Reduce the expression 

X ^ ci cos kt + ctonkt 

to the form 

X ^ Acoaikt + e). 

Find A and in terms of Ci and Cs. 

2. Reduce the expression 

X =8 er^i [ci cos /3< + Cj sin /3fl 
to the form (33b). 
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8. The force exerted by a spring is proportional to the amount the 
spring is stretched and i 3 2 lbs. when the spring is stretched i in. A 
mass of 6 lbs. is suspenc ed by the spring. If the mass is depressed 
slightly and released, find the time of vibration. 

4. A disk of radius a and' mass M is supported in a horizontal plane 
by a vertical wire attached to its center. When the disk is rotated 
about the vertical axis through the angle d the torque exerted by the 
wire is 

When the disk is rotated through a small angle and released it vibrates 
n times per second. Find the constant k. 

6. A cylindrical spar buoy stai^ds vertically in the water. Its 
diameter is 12 in. and its mass 200 lbs. Find the time of vibration 
when it is depressed sUghtly and released. 

6. Two equal weights are hanging a^ the end of an elastic string. 
One falls off*. Find the motion of the i^'emaining weight. 

7. One end of an elastic string of natur il length a is attached to a 
fixed point on a horizontal table and to ^he other end a mass m is 
attached. If m is drawn aside until the strL \g is elongated the amount 
a and released, find the time of one comple1;e vibration if a force of 
P lbs. elongates the string 1 ft. 

8. A 10-lb. body is observed to vibrate 90 times per minute and the 
oscillation damps to } amplitude in 15 sees. Find the differential 
equation of motion. 

9. A 10-lb. mass is acted upon by a restitutive force proportional to 
the displacement and equal to 2 lbs. when the displacement is 1 ft. 
Find the period if the vibration damps to f^r amphi'ude after 3 com- 
plete vibrations. 

10. A rod of length 2 lis supported in a horizontal position by two 
vertical strings of length I attached to its ends. If the rod is turned 
through a small angle about the vertical line through its center and 
released, find the time of vibration. 

11. A rod of length 2 1 rests in a horizontal position on a cylinder of 
radius a. If it is tipped sUghtly, find the time of vibration. 

^12. If C is the capacity of a condenser the electromotive force across 
the condenser is 

idt. 



if' 



If a constant electromotive force E is impressed on a circuit containing 
a resistance Ry inductance L, and capacity C, find the current as a 

function of the time if i == and I i dt « when t ™ 0. (See Prob. 

3, page 31.) 



/ 
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w 

18. A wire of length I and mass m is f astenedi by one end and rotated 
with angular velocity w in a horizontal plane. If under a force F 
it stretches the amount IdFy find the amount it is stretched by centrif- 
ugal force. 

14. A horizontal tube rotates about a vertical axis with angular 
velocity w. A ball inside the tube and shading without friction starts 
from the center with velocity t;©. Find its motion. 

16. A particle of mass m is attracted, toward each of two centers of 
force with a force equal to k times tbje distance. If the distance be- 
tween the two centers is 2 & and thfj particle starts from rest at the 
point on the joining line at distance, c from the center, find its motion. 

16. In the preceding problem fijid the motion if the particle starts 
with velocity v^ perpendicular to the line joining the centers of force, 

17. When a current i flows tjhrough a resistance R the drop in po- 
tential is %R. When a potentifiil E is impressed on a region of leakance 
Gf the current that leaks away is EG. If the resistance of a cable is r 
per unit length and its leakjince g per unit length, what current and 
potential must a line of length I receive at one end if the current and 
potential at the other end are to be /o and E^, 

18. A helical spring of natural length 21 and negligible mass is 
hung up by one end. Two equal weights of mass m are attached, one 
at the middle, the other at the lower end, and let fall. Find the motion 
if when hanging at rest the weights stretch the spring the amount c, 

19. A weight of 4 lbs. is suspended by a spring. The weight extends 
the spring 1 inch. If the upper end of the spring is given the harmonic 
motion 

y = sin {Vl2g t) 

in the vertical line, find the motion of the suspended weight.- 

20. Two weights of 4 lbs. are connected by the spring of the pre- 
ceding problem, one of the weights being held fixed and the other 
suspended by the spring. If the upper weight is released, find the 
motion of the two weights. 

21. According to Newton's law the gravitational attraction between 
two masses mi, mx at distance d is 

Find the motion of the two masses assmning that their center of cpeav- 
ity remains fixed. 



ANSWERS 



1. tan^ — cotV " c. 

2. 1 + 2/* = c (1 •- a?). 

3. 3^+3? -^ = e. 

4. y + a = c sin X. 



Pages 3, 4 

5. X 



= ln 



6. a; = 



ev — 1 

c» • 

V 



Pages 14-18 

2 - vr^^ 



1. ^ = car. 

2. y - V4 -x» + 2 ln( 

y = 6.662/0 4. i2 = /2a c-.l78» 

(0 = 100 e-l-l* revolutions per minute. 



X 



) 



3. 

5. 

6. 

7. 

9. 
10. 
11. 
12. 
13. 
14. 
15. 
16. 
17. 



Ihr. 

Xo = 56.5. 7.84 hrs. 

Parabola. 

Paraboloid of revolution. 

Hyperboloid of revolution. 

10 years. 

About 11} years. 

iJP + hW) Td, 
p = 14.7 e-.oooo4» 
About 17} miles. 



19. T. _ ^ = f r. _ !£?! W 

\ 9 / 



22. 9.8 min. 

23. 10.2 min. 

24. 4 min. 

25. 18.4 

26. 4.6 min. 

27. 6.6 days. 

28. 19.5 lbs. 

29. 54.7 lbs. 

30. 99.5 %. 

31. 1500 cu. ft. 

32. .24%. 

33. .124%. 

34. T = ix. 864,000 cal. 



35. 
36. 

^7. 

38. 

1. 
2. 
3. 
4. 



r « 592 - 187.6 In r. 1,731,000 cal. 
.0005 39. 2.9 days. 

1 ^ 40. 15.5 lbs. 

1024' 41. 50 min. 

2.»4bs 



x^ + a^ — xy^ 
x^ ^ Sxy + c. 
aj* + 3 a^l/" = c. 
a^ — 2/* = ca;. 



Pages 29, 80 

— 2/* = c. 5. 2xy -\-^ - c. 

6. 

7. 2/ (1 + ex) = 1 
71 



In 2/ = c. 

2/ 
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8. » - j^ + ce«. 16. -i = X + i + e«*. 

9. «=-^«.'-l). ''■ *'=''(*^+'')- 

10. ,=^8«. !«• ln(^+2/')=2tan-tf + c. 

X 19. J/* - 1 = c (1 + x)«e-2x. 

12. y+a = csma;. ^' ^"2^ "27; 

13. X* + 4 1/ (a;* - 1)1 = c 21. 1/ sin X + i 2^ - i X* = c. 

14. X = !/• (c - e-v). 22. y = J x [e*+c _ c-«-^]. 

15. ev + In X = c. 

23. (2 y - x) Vl + x» = c - In (x + Vx« + 1). 

24. y* = 4 xy + c. 25. 2/* = } sin X + c esc* x. 

26. x+22/+ln(x+2/ -2)'=c. 

27. y» = ce* - X - 1. 29. y = ce^f or 2/ = c + i x*. 

. C . 1 • 30. e2a: - 2 ea^+1f = c. 

28. 2/=2^-2V 

31. X + 2 2/ + 3 In (2 X + 3 2/ - 7) = c. 

32. (x +2/)* = cx«e-^. 

33. X = i e« + cie-', j/ = i c« - cie-* + Ct. 

34. X = cie-' + C2e' + 1, 2/ = cie-' — C2e' + 1. 

35. y = ciel' + cje**, x « 4 ciel< + 3 C26i«. 

36. X = cie-« + dfi-^, y = Cie-* + 3 c^-^ + cos «. 

Pages 31, 32 

2. 14.9 lbs. 

E r --H 

3. i = v>« . y « a I i2 sin co< — I/O) cos ca< + Zxoe /' I* 



/? + 



4. ii = 



p-^^ r(/2i22 + LL2f^) sin co< 



/2* + L*«2 



_A -1 
+ oi{L%R — lUL) (cos orf — e i ) , 

where -R = i2i + /22, L =^ Li+L%. MS^M!' 

Wtt yM 5. rp = c (eJ*^'^ — 1), c, A; being constants.'^ 1 1^ 6 •• 
^6. x = i(l -2-'),2/ = i(l -2-0. 

7. 2/ = 4.5 (1 - e-.7324«), a? = 10 - 2/. 

8. B = .214, c - .249. g. (i+?)sind = c. 

Page 36 

1. 2/ = ^ In X + i x* + cix + C2. 4. y = ciex + C2€-*. 

2. 2/ = Ci In X + C2 — i x^. 5. y = ci cos (fcx + Cj). 

3. 2/ = i (^ + «)' + cix + cj. 



ANSWERS 73 



6. ci<+C2 = V«(ci8+2A;«)+-^8in-i(^^^^^ 

7. efl«» = ^ [ef^f-^i 4- e-o^'-^a). 8. y ^ i cix^ - J^hi x + d, 

2 Z Ci 

9. 2/ = a; 4- cixea? + d. 11. 2/ = J On a;)* + ci In a; + Cj. 

10. y = i (eCiaj+«2 4- e-ci^^-cs). 12. ciS* = a« + (ci< + cs)«. 



Pages 61-64 

1. Ely = 1 TT (te* — i X*), the origin being the middle point 

2. Ely = i TT (i a;* — ix^), the fixed end being the origiA. 

3. Ely = -^w;(6Pa;*-42a:»4- a;*). 

4. J^/t/ =^w;[3Pa:»-2te«- a;^i; the length being 2 Z. 

5. The same as Prob. 4 if the length is I, 

6. EIy=:^l2Px'-n 

^- ^"V^ ,' 10. t/ = c(et* + e-to) 

11. y = - f ea + e"« j, if the axes are properly chosen. 

12. 65.5 lbs. 4/ n^ 

13. Parabola. ^^' ^'\ i 



mg 4- kv^, 
1^- F' 23. « = y ^hi (6 4- v^). 



15. 77.6 ft. 2.15 sec. 



g 



17. About 7 miles per second. 24. i = -^ In (9 + 4 Vs). 

18. 116 hours. \^ 

19. About 42i mm. 25. < = -^hi (17 4- 4 vTs). 

20. 8 = i flf (sin a — M cos a) ^. 

26. Maximum height 48,500 ft., range 36.7 mi. 

27. a? = T (1 - e m'), 

y = ^ (fcvo sm a 4- w^) (1 — e m J — ~- • 

28. ^+-^. = 1. ' 29. ^-^-1. 

30. Helix. 
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32. 'nieempse,l = (i-^.)co8^ + 



33. r = 



2 To 34. r s= a cos 0. 



1 + coae 
m / . 2K\ /. I'J^N 2^ 



36. 77 (y<* sin 0. 
14 

2 w 

37. It slides until i = =— ^ and then rolls. 

4?iui Oft 2mji 

40. The distance M moves in t seconds is 



2m+4Af 



Pages G6, 66 

1. 2/ = ci + cje*. 4. 2/ == ci cos X + C2 sin a?. 

2. y = cieaf + cse*«f. 6. y = Ci + Cjc-a; + c^e**. 

3. t/ = (ci + cjx)e**. 

6. y = Cie^f + C2e-* + cs sin a; + C4 cos x. 

7. y ^ &^ [ci cos X + C2 sin xj. 

8. 2/ = e-Ja; [ci cos (i V3xy*4^C2 sm (iv^ x),\ 

9. 2/ = (ci 4- Cax + CzPi^)e-x, 

10. y == 2 — a; + ci cos x + C2 sin x. 

11. y = cie** + C2fir-^ - i a:* - *. 

12. 2/ = cie^ — i (sin x + cos x). 

13. 2/ = Ci + cjeaf - i X* - X. 

4 1 

14. 2/ = cie-^ + C2e-»^ + *^"9+i5^- 

15. y = cie«« 4- C3fi-ax 4- __ aje<«. 

16. 1/ = ( Ci + s— ) sin ax 4- C2 cos ax, 

, / X V3 . . X Vs\ ^ c 

17. 2/ = cie» 4- e^^ ( ca cos— ^ h dsm—^J — ^ — f^ 

18. y = Cie^f 4- Cac*^ — J e** sin x. 
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19. y = Ci6»* + cje-«« +5= 6"* (6 sin x - cos x). 

20. y = ci + CiX + cax» + Cie^ + j^gg (4 cos 4 a; - sin 4 x). 

21. y--(ci+w + ia^)e-^ + i c*. 

22. y = ci^ + C4e-< — 1, x = cie< — cjc-' — 1. 

23. y = ci cos < + Ci sin «, x = (ci + cj) cos < + (c* — <?i) sin t 

24. X « Cie-< + cje-*', y = cie-< + 3 oje-"' + cx)s <. 

25. X = cie* + C2e~' + cj cos < + C4 sin <, 
y = ci€^ + C4«^ — Ci cos t — CiOnt, 

26. X = ci + C4< + C3«* - J ^ + e« 

y = C4 - (ci +2c») < - J (C4 - 1) <• - Jci^ +4<* - e*. 



3 



Pages 68-70 

• 2V FS' 



2V 3 ^ 5* 2.24 sec. 

6. x = acos(y-n, where a is the amount the string is stretched 

by one weight and x is measured from the point where it would hang 
in equilibrium. 

7. (4+2x)\/|- 

8. ^ + .0924^ + 88.83x = 0. 

9. 2.49 sec. 

2irl 



i 



10. 2.\/±. 11. 



3j' VSop 

12. i- ^^^ ^§i^„, ViLC-C^E^ 



13. / (— 1 ) , where ^ = w y/^^^ 



14. f = i5L (e-« - e-f). 
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15. X = c cos ^V — , the origin being midway between the centers. 

16. An ellipse. 

17. i - j(/. + E,\f«>f^^' + i(/. - B. V^) 
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